A study of transposable element-associated structural variations (TASVs) using a de novo-assembled Korean genome.

Advances in next-generation sequencing (NGS) technology have made personal genome sequencing possible, and indeed, many individual human genomes have now been sequenced. Comparisons of these individual genomes have revealed substantial genomic differences between human populations as well as between individuals from closely related ethnic groups. Transposable elements (TEs) are known to be one of the major sources of these variations and act through various mechanisms, including de novo insertion, insertion-mediated deletion, and TE-TE recombination-mediated deletion. In this study, we carried out de novo whole-genome sequencing of one Korean individual (KPGP9) via multiple insert-size libraries. The de novo whole-genome assembly resulted in 31,305 scaffolds with a scaffold N50 size of 13.23 Mb. Furthermore, through computational data analysis and experimental verification, we revealed that 182 TE-associated structural variation (TASV) insertions and 89 TASV deletions contributed 64,232 bp in sequence gain and 82,772 bp in sequence loss, respectively, in the KPGP9 genome relative to the hg19 reference genome. We also verified structural differences associated with TASVs by comparative analysis with TASVs in recent genomes (AK1 and TCGA genomes) and reported their details. Here, we constructed a new Korean de novo whole-genome assembly and provide the first study, to our knowledge, focused on the identification of TASVs in an individual Korean genome. Our findings again highlight the role of TEs as a major driver of structural variations in human individual genomes.

Epigenetic analysis in rheumatoid arthritis synoviocytes.

Rheumatoid arthritis (RA) is a complex chronic systematic disease with progressive destruction of the joints by invasive synoviocytes. To characterize the key regulators involved in the development of RA, we obtained multilayer epigenomics data including DNA methylation by whole-genome bisulfite sequencing, miRNA profiles, genetic variations by whole-exome sequencing, and mRNA profiles from synoviocytes of RA and osteoarthritis (OA) patients. The overall DNA methylation patterns were not much different between RA and OA, but 523 low-methylated regions (LMRs) were specific to RA. The LMRs were preferentially localized at the 5' introns and overlapped with transcription factor binding motifs for GLI1, RUNX2, and TFAP2A/C. Single base-scale differentially methylated CpGs were linked with several networks related to wound response, tissue development, collagen fibril organization, and the TGF-ß receptor signaling pathway. Further, the DNA methylation of 201 CpGs was significantly correlated with 27 expressed miRNA genes. Our interpretation of epigenomic data of the synoviocytes from RA and OA patients is an informative resource to further investigate regulatory elements and biomarkers responsible for the pathophysiology of RA and OA.

Differential DNA methylation of MSI2 and its correlation with diabetic traits.

Differential DNA methylation with hyperglycemia is significantly associated with Type 2 Diabetes (T2D). Longtime extended exposure to high blood glucose levels can affect the epigenetic signatures in all organs. However, the relevance of the differential DNA methylation changes with hyperglycemia in blood with pancreatic islets remains unclear. We investigated differential DNA methylation in relation to glucose homeostasis based on the Oral Glucose Tolerance Test (OGTT) in a population-based cohort. We found a total of 382 differential methylation sites from blood DNA in hyperglycemia and type 2 diabetes subgroups using a longitudinal and cross-sectional approach. Among them, three CpG sites were overlapped; they were mapped to the MSI2 and CXXC4 genes. In a DNA methylation replication study done by pyrosequencing (n = 440), the CpG site of MSI2 were shown to have strong associations with the T2D group (p value = 2.20E-16). The differential methylation of MSI2 at chr17:55484635 was associated with diabetes-related traits, in particular with insulin sensitivity (QUICKI, p value = 2.20E-16) and resistance (HOMA-IR, p value = 1.177E-07). In human pancreatic islets, at the single-base resolution (using whole-genome bisulfite sequencing), the 292 CpG sites in the ±5kb at chr17:55484635 were found to be significantly hypo-methylated in donors with T2D (average decrease = 13.91%, 95% confidence interval (CI) = 4.18~ 17.06) as compared to controls, and methylation patterns differed by sex (-9.57%, CI = -16.76~ -6.89) and age (0.12%, CI = -11.17~ 3.77). Differential methylation of the MSI2 gene (chr17:55484635) in blood and islet cells is strongly related to hyperglycemia. Our findings suggest that epigenetic perturbation on the target site of MSI2 gene in circulating blood and pancreatic islets should represent or affect hyperglycemia.

Mitochondrial DNA copy number augments performance of A1C and oral glucose tolerance testing in the prediction of type 2 diabetes.

Here, we tested the performance of the mitochondrial DNA copy number (mtDNA-CN) in predicting future type 2 diabetes (n = 1108). We used the baseline clinical data (age, sex, body mass index, waist-to-hip ratio, systolic and diastolic blood pressure) and the mtDNA-CN, hemoglobin A1c (A1C) levels and results of oral glucose tolerance test (OGTT) including fasting plasma glucose, 1-hour glucose, and 2-hour glucose levels, to predict future diabetes. We built a prediction model using the baseline data and the diabetes status at biannual follow-up of 8 years. The mean area under curve (AUC) for all follow-ups of the full model including all variables was 0.92 ± 0.04 (mean ± standard deviation), while that of the model excluding the mtDNA-CN was 0.90 ± 0.03. The sensitivity of the f4ull model was much greater than that of the model not including mtDNA-CN: the mean sensitivities of the model with and without mtDNA-CN were 0.60 ± 0.06 and 0.53 ± 0.04, respectively. We found that the mtDNA-CN of peripheral leukocytes is a biomarker that augments the predictive power for future diabetes of A1C and OGTT. We believe that these results could provide invaluable information for developing strategies for the management of diabetes.

Genome-Wide Association Study Meta-Analysis of Long-Term Average Blood Pressure in East Asians.

BACKGROUND: Genome-wide single marker and gene-based meta-analyses of long-term average (LTA) blood pressure (BP) phenotypes may reveal novel findings for BP. METHODS AND RESULTS: We conducted genome-wide analysis among 18 422 East Asian participants (stage 1) followed by replication study of ≤46 629 participants of European ancestry (stage 2). Significant single-nucleotide polymorphisms and genes were determined by a P<5.0×10-8 and 2.5×10-6, respectively, in joint analyses of stage-1 and stage-2 data. We identified 1 novel ARL3 variant, rs4919669 at 10q24.32, influencing LTA systolic BP (stage-1 P=5.03×10-8, stage-2 P=8.64×10-3, joint P=2.63×10-8) and mean arterial pressure (stage-1 P=3.59×10-9, stage-2 P=2.35×10-2, joint P=2.64×10-8). Three previously reported BP loci (WBP1L, NT5C2, and ATP2B1) were also identified for all BP phenotypes. Gene-based analysis provided the first robust evidence for association of KCNJ11 with LTA systolic BP (stage-1 P=8.55×10-6, stage-2 P=1.62×10-5, joint P=3.28×10-9) and mean arterial pressure (stage-1 P=9.19×10-7, stage-2 P=9.69×10-5, joint P=2.15×10-9) phenotypes. Fourteen genes (TMEM180, ACTR1A, SUFU, ARL3, SFXN2, WBP1L, CYP17A1, C10orf32, C10orf32-ASMT, AS3MT, CNNM2, and NT5C2 at 10q24.32; ATP2B1 at 12q21.33; and NCR3LG1 at 11p15.1) implicated by previous genome-wide association study meta-analyses were also identified. Among the loci identified by the previous genome-wide association study meta-analysis of LTA BP, we transethnically replicated associations of the KCNK3 marker rs1275988 at 2p23.3 with LTA systolic BP and mean arterial pressure phenotypes (P=1.27×10-4 and 3.30×10-4, respectively). CONCLUSIONS: We identified 1 novel variant and 1 novel gene and present the first direct evidence of relevance of the KCNK3 locus for LTA BP among East Asians.

Genome-based exome sequencing analysis identifies GYG1, DIS3L and DDRGK1 are associated with myocardial infarction in Koreans.

Myocardial infarction (MI) is a complex disease caused by combination of genetic and environmental factors. Although genome-wide association studies (GWAS) identified more than 46 risk loci which are associated with coronary artery disease and MI, most of the genetic variability inMI still remains undefined. Here, we screened the susceptibility loci for MI using exome sequencing and validated candidate variants in replication sets. We identified that three genes (GYG1, DIS3L and DDRGK1) were associated with MI at the discovery and replication stages. Further research will be required to determine the functional association of these genes with MI risk, and these associations have to be confirmed in other ethnic populations.

Inorganic Phosphorus and Potassium Are Putative Indicators of Delayed Separation of Whole Blood.

OBJECTIVES: The delayed separation of whole blood can influence the concentrations of circulating blood components, including metabolites and cytokines. The aim of this study was to determine whether clinical-biochemistry analytes can be used to assess the delayed separation of whole blood. METHODS: We investigated the plasma and serum concentrations of five clinical-biochemistry analytes and free hemoglobin when the centrifugation of whole blood stored at 4°C or room temperature was delayed for 4 hours, 6 hours, 24 hours, or 48 hours, and compared the values with those of matched samples that had been centrifuged within 2 hours after whole-blood collection. RESULTS: The inorganic phosphorus (IP) levels in the plasma and serum samples were elevated ≥ 1.5-fold when whole-blood centrifugation was delayed at room temperature for 48 hours. Furthermore, the IP levels in the plasma samples showed excellent assessment accuracy [area under the receiver-operating-characteristic curve (AUC) > 0.9] after a 48-hour delay in whole-blood separation, and high sensitivity (100%) and specificity (95%) at an optimal cutoff point. The IP levels in the serum samples also exhibited good assessment accuracy (AUC > 0.8), and high sensitivity (81%) and specificity (100%). The potassium (K(+)) levels were elevated 1.4-fold in the serum samples following a 48-hour delay in whole-blood separation. The K(+) levels showed excellent assessment accuracy (AUC > 0.9) following a 48-hour delay in whole-blood separation, and high sensitivity (95%) and specificity (91%) at an optimal cutoff point. CONCLUSION: Our study showed that the IP and K(+) levels in the plasma or serum samples could be considered as putative indicators to determine whether whole-blood separation had been delayed for extended periods.

Genome-wide association studies in East Asians identify new loci for waist-hip ratio and waist circumference.

Sixty genetic loci associated with abdominal obesity, measured by waist circumference (WC) and waist-hip ratio (WHR), have been previously identified, primarily from studies conducted in European-ancestry populations. We conducted a meta-analysis of associations of abdominal obesity with approximately 2.5 million single nucleotide polymorphisms (SNPs) among 53,052 (for WC) and 48,312 (for WHR) individuals of Asian descent, and replicated 33 selected SNPs among 3,762 to 17,110 additional individuals. We identified four novel loci near the EFEMP1, ADAMTSL3 , CNPY2, and GNAS genes that were associated with WC after adjustment for body mass index (BMI); two loci near the NID2 and HLA-DRB5 genes associated with WHR after adjustment for BMI, and three loci near the CEP120, TSC22D2, and SLC22A2 genes associated with WC without adjustment for BMI. Functional enrichment analyses revealed enrichment of corticotropin-releasing hormone signaling, GNRH signaling, and/or CDK5 signaling pathways for those newly-identified loci. Our study provides additional insight on genetic contribution to abdominal obesity.

High-density genotyping of immune-related loci identifies new SLE risk variants in individuals with Asian ancestry.

Systemic lupus erythematosus (SLE) has a strong but incompletely understood genetic architecture. We conducted an association study with replication in 4,478 SLE cases and 12,656 controls from six East Asian cohorts to identify new SLE susceptibility loci and better localize known loci. We identified ten new loci and confirmed 20 known loci with genome-wide significance. Among the new loci, the most significant locus was GTF2IRD1-GTF2I at 7q11.23 (rs73366469, Pmeta = 3.75 × 10(-117), odds ratio (OR) = 2.38), followed by DEF6, IL12B, TCF7, TERT, CD226, PCNXL3, RASGRP1, SYNGR1 and SIGLEC6. We identified the most likely functional variants at each locus by analyzing epigenetic marks and gene expression data. Ten candidate variants are known to alter gene expression in cis or in trans. Enrichment analysis highlights the importance of these loci in B cell and T cell biology. The new loci, together with previously known loci, increase the explained heritability of SLE to 24%. The new loci share functional and ontological characteristics with previously reported loci and are possible drug targets for SLE therapeutics.

Impact of Whole-Blood Processing Conditions on Plasma and Serum Concentrations of Cytokines.

Pre-analytical variations in plasma and serum samples can occur because of variability in whole-blood processing procedures. The aim of this study was to determine the impact of delayed separation of whole blood on the plasma and serum concentrations of cytokines. The concentrations of 16 cytokines were measured in plasma and serum samples when the centrifugation of whole blood at room temperature was delayed for 4, 6, 24, or 48 h, and the values were compared with those observed after separation within 2 h of whole-blood collection. Receiver operating characteristic (ROC) curve analysis was also performed for cytokines to determine whether cytokine levels in plasma and serum samples can be used to assess delayed separation of whole blood. Plasma concentrations of interleukin (IL)-1ß, granulocyte-macrophage colony-stimulating factor (GM-CSF), and soluble CD40 ligand (sCD40L) and serum concentrations of IL-1ß, IL-6, IL-8, macrophage inflammatory protein-1α (MIP-1α), and MIP-1ß increased significantly (>2-fold) when separation was delayed at room temperature for 24 h. The concentrations of 6 of these cytokines (all except serum IL-1ß and IL-6) demonstrated high diagnostic performance (area under the ROC curve >0.8) for delayed separation of whole blood. Furthermore, these cytokine concentrations typically exhibited high sensitivity and specificity at each optimal cutoff point. Conversely, IL-17A was stable in both plasma and serum samples, even when whole-blood centrifugation was delayed at room temperature for 48 h. This study shows that certain cytokines (IL-1ß, GM-CSF, sCD40L, IL-8, MIP-1α, and MIP-1ß) could be used for assessing the quality of plasma or serum samples.

Meta analysis identifies a novel susceptibility locus associated with heel bone strength in the Korean population.

INTRODUCTION: Calcaneal quantitative ultrasound has been recognized as a non-invasive method for evaluation of bone strength and prediction of osteoporotic fracture. METHODS: To extend a thorough genetic catalog for osteoporotic bone properties, we performed a genome-wide association study (rural cohort I, n=1895) of speed of sound (SOS) using the 1000 genome-based imputation in the discovery stage and then carried out in silico lookups (rural cohort II and III, n=2,967) and de novo genotyping (rural cohort IV, n=4,296) in the replication stage. RESULTS: In the combined meta-analysis (n=9,158), we identified a novel variant associated with SOS (rs2445771 in the GLDN gene, P=2.27×10(-9)) reaching genome-wide significance in the Korean population. We further demonstrated that allele-specific regulatory modifications found to be associated with functional enrichments by ENCODE annotations. CONCLUSION: Our findings could provide additional insights into understanding of genetic and epigenetic regulations on bone metabolism.

An epigenomic signature of postprandial hyperglycemia in peripheral blood leukocytes.

Postprandial hyperglycemia is known to be one of the earliest signs of abnormal glucose homeostasis associated with type 2 diabetes. This study aimed to assess clinical significance of a 1-h postprandial glucose level for the development of diabetes, and identify epigenetic biomarkers of postprandial hyperglycemia. We analyzed clinical data from the oral glucose tolerance tests for healthy subjects (n=4502). The ratio (Glu60/Glu0) of 1-h glucose levels to fasting glucose levels was significantly associated with an insulin sensitive index (QUICKI, quantitative insulin sensitivity check index) (ß=0.055, P=1.25E-04) as well as a risk of future pre-diabetic and diabetic conversion. Next, DNA methylation profile analyses of 24 matched pairs of the high and low Glu60/Glu0 ratio subjects showed that specific DNA methylation levels in the promoter region of an olfactory receptor gene (olfactory receptor gene family10 member A4, OR10A4) were associated with the Glu60/Glu0 ratios (ß=0.337, P=0.03). Moreover, acute oral glucose challenges decreased the DNA methylation levels of OR10A4 but not the global DNA methylation in peripheral leukocytes of healthy subjects (n=7), indicating that OR10A4 is a specific epigenomic target of postprandial hyperglycemia. This work suggests possible relevance of olfactory receptor genes to an earlier molecular biomarker of peripheral hyperglycemia and diabetic conversion.

Contribution of GABRG2 Polymorphisms to Risk of Epilepsy and Febrile Seizure: a Multicenter Cohort Study and Meta-analysis.

Gamma-aminobutyric acid receptor (GABA-A) is the most common receptor of fast synaptic inhibition in the human brain. Gamma protein encoded by the GABRG2 gene is one of the subunits of the GABA-A receptor, which plays an essential role in the function of this receptor. Several studies have identified various febrile seizure (FS) and epilepsy risk variants of GABRG2 gene in different populations, but some others did not support these results. The aim of this case-control study is to investigate whether GABRG2 polymorphisms contribute to susceptibility for FS and epilepsy in pooled data of three cohorts, from Malaysia (composed of Malay, Chinese, and Indian), Hong Kong, and Korea. Furthermore, the pooled dataset of these cohorts with previous reports were meta-analyzed for determining the risk effect size of the rs211037 polymorphism on FS and symptomatic epilepsy (SE). The rs211037, rs210987, rs440218, rs2422106, rs211014, and rs401750 polymorphisms were genotyped in the 6442 subjects (1729 epilepsy and 4713 controls). Results of the case-control study showed associations between rs211037 and the risk of SE in the pooled data from all cohorts (T vs. C, p = 3 × 10(-5), and TT vs. CC, p = 2 × 10(-5)) and the risk of partial seizure in the combined data of Malaysia and Hong Kong (both T vs. C and TT vs. CC, p = 2 × 10(-6)). The rs211037-rs210987 and rs2422106-rs211014-rs401750 haplotypes were also associated with susceptibility to SE in Chinese. Meta-analysis of all Asians identified association between rs211037 and FS and SE (T vs. C, p = 4 × 10(-4), and p = 4 × 10(-3), respectively). In conclusion, rs211037 alone may be a risk factor for FS, partial seizure, and SE, and in linkage disequilibrium with rs210987 can contribute to FS and SE in Asians, particularly in Chinese.

Identification of a Systemic Lupus Erythematosus Risk Locus Spanning ATG16L2, FCHSD2, and P2RY2 in Koreans.

OBJECTIVE: Systemic lupus erythematosus (SLE) is a chronic autoimmune disorder whose etiology is incompletely understood, but likely involves environmental triggers in genetically susceptible individuals. Using an unbiased genome-wide association (GWA) scan and replication analysis, we sought to identify the genetic loci associated with SLE in a Korean population. METHODS: A total of 1,174 SLE cases and 4,246 population controls from Korea were genotyped and analyzed with a GWA scan to identify single-nucleotide polymorphisms (SNPs) significantly associated with SLE, after strict quality control measures were applied. For select variants, replication of SLE risk loci was tested in an independent data set of 1,416 SLE cases and 1,145 population controls from Korea and China. RESULTS: Eleven regions outside the HLA exceeded the genome-wide significance level (P = 5 × 10(-8) ). A novel SNP-SLE association was identified between FCHSD2 and P2RY2, peaking at rs11235667 (P = 1.03 × 10(-8) , odds ratio [OR] 0.59) on a 33-kb haplotype upstream of ATG16L2. In the independent replication data set, the SNP rs11235667 continued to show a significant association with SLE (replication meta-analysis P = 0.001, overall meta-analysis P = 6.67 × 10(-11) ; OR 0.63). Within the HLA region, the SNP-SLE association peaked in the class II region at rs116727542, with multiple independent effects observed in this region. Classic HLA allele imputation analysis identified HLA-DRB1*1501 and HLA-DQB1*0602, each highly correlated with one another, as most strongly associated with SLE. Ten previously established SLE risk loci were replicated: STAT1-STAT4, TNFSF4, TNFAIP3, IKZF1, HIP1, IRF5, BLK, WDFY4, ETS1, and IRAK1-MECP2. Of these loci, previously unreported, independent second risk effects of SNPs in TNFAIP3 and TNFSF4, as well as differences in the association with a putative causal variant in the WDFY4 region, were identified. CONCLUSION: Further studies are needed to identify true SLE risk effects in other loci suggestive of a significant association, and to identify the causal variants in the regions of ATG16L2, FCHSD2, and P2RY2.

MicroRNA-650 in a copy number-variable region regulates the production of interleukin 6 in human osteosarcoma cells.

Copy number variation is a well-known genetic variation. microRNAs (miRNAs/miRs) are non-coding RNAs that mediate gene expression by regulating target mRNAs. In the present study, copy number deletions encompassing miRNA coding regions were investigated to determine the association between the deletion of miRNA and its phenotypic effects. A total of 38 human miRNAs in copy number variants were identified and miR-650, which is functional in the human osteosarcoma MG-63 cell line, was selected. Overexpression of miR-650 decreased the expression of inhibitor of growth family member 4 (ING4) in the MG-63 cells and increased interleukin (IL)6 transcription, as well as IL6 secretion in IL1B-stimulated cells. Furthermore, miR-650 downregulated the amount of nuclear factor of κ light polypeptide gene enhancer in B cells inhibitor α and increased the transcriptional activity of nuclear factor (NF)κB. Downregulation of ING4 also increased the production of IL6, similar to miR-650 overexpression. Taken together, these data indicate that miR-650 plays a significant role in the production of IL6 by regulating ING4 expression and NFκB signaling in IL1B-stimulated MG-63 osteosarcoma cells.

On the Estimation of Heritability with Family-Based and Population-Based Samples.

For a family-based sample, the phenotypic variance-covariance matrix can be parameterized to include the variance of a polygenic effect that has then been estimated using a variance component analysis. However, with the advent of large-scale genomic data, the genetic relationship matrix (GRM) can be estimated and can be utilized to parameterize the variance of a polygenic effect for population-based samples. Therefore narrow sense heritability, which is both population and trait specific, can be estimated with both population- and family-based samples. In this study we estimate heritability from both family-based and population-based samples, collected in Korea, and the heritability estimates from the pooled samples were, for height, 0.60; body mass index (BMI), 0.32; log-transformed triglycerides (log TG), 0.24; total cholesterol (TCHL), 0.30; high-density lipoprotein (HDL), 0.38; low-density lipoprotein (LDL), 0.29; systolic blood pressure (SBP), 0.23; and diastolic blood pressure (DBP), 0.24. Furthermore, we found differences in how heritability is estimated--in particular the amount of variance attributable to common environment in twins can be substantial--which indicates heritability estimates should be interpreted with caution.

Global implementation of genomic medicine: We are not alone.

Around the world, innovative genomic-medicine programs capitalize on singular capabilities arising from local health care systems, cultural or political milieus, and unusual selected risk alleles or disease burdens. Such individual efforts might benefit from the sharing of approaches and lessons learned in other locales. The U.S. National Human Genome Research Institute and the National Academy of Medicine recently brought together 25 of these groups to compare projects, to examine the current state of implementation and desired near-term capabilities, and to identify opportunities for collaboration that promote the responsible practice of genomic medicine. Efforts to coalesce these groups around concrete but compelling signature projects should accelerate the responsible implementation of genomic medicine in efforts to improve clinical care worldwide.

Prediction of breast cancer survival using clinical and genetic markers by tumor subtypes.

PURPOSE: To identify the genetic variants associated with breast cancer survival, a genome-wide association study (GWAS) was conducted of Korean breast cancer patients. METHODS: From the Seoul Breast Cancer Study (SEBCS), 3,226 patients with breast cancer (1,732 in the discovery and 1,494 in the replication set) were included in a two-stage GWAS on disease-free survival (DFS) by tumor subtypes based on hormone receptor (HR) and human epidermal growth factor receptor 2 (HER2). The associations of the re-classified combined prognostic markers through recursive partitioning analysis (RPA) of DFS for breast cancer were assessed with the Cox proportional hazard model. The prognostic predictive values of the clinical and genetic models were evaluated by Harrell's C. RESULTS: In the two-stage GWAS stratified by tumor subtypes, rs166870 and rs10825036 were consistently associated with DFS in the HR+ HER2- and HR- HER2- breast cancer subtypes, respectively (Prs166870 = 2.88 × 10(-7) and Prs10825036 = 3.54 × 10(-7) in the combined set). When patients were classified by the RPA in each subtype, genetic factors contributed significantly to differentiating the high risk group associated with DFS inbreast cancer, specifically the HR+ HER2- (P discovery=1.18 × 10(-8) and P replication = 2.08 × 10(-5)) and HR- HRE2- subtypes (P discovery = 2.35 × 10(-4) and P replication = 2.60 × 10(-2)). The inclusion of the SNPs tended to improve the performance of the prognostic models consisting of age, TNM stage and tumor subtypes based on ER, PR, and HER2 status. CONCLUSION: Combined prognostic markers that include clinical and genetic factors by tumor subtypes could improve the prediction of survival in breast cancer.